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Abstract 

We investigate the phenomenon of decay of a supercurrent in a superconducting 
thin film in the absence of an applied magnetic field. The resulting zero-temperature 
resistance derives from two equally possible mechanisms: 1) quantum tunneling of 
vortices from the edges of the sample; and 2) homogeneous quantum nucleation of 
vortex-antivortex pairs in the bulk of the sample, arising from the instability of the 
Magnus field's "vacuum" . We study both situations in the case where quantum dissi- 
pation dominates over the inertia of the vortices. We find that the vortex tunneling and 
nucleation rates have a very rapid dependence on the current density driven through 
the sample. Accordingly, whilst normally the superconductor is essentially resistance- 
free, for the high current densities that can be reached in high-T c films a measurable 
resistance might develop. We show that edge-tunneling appears favoured, but the 
presence of pinning centres and of thermal fluctuations leads to an enhancement of 
the nucleation rates. In the case where a periodic pinning potential is artificially in- 
troduced in the sample, we show that current-oscillations will develop indicating an 
effect specific to the nucleation mechanism where the vortex pair-production rate, thus 
the resistance, becomes sensitive to the corrugation of the pinning substrate. In all 
situations, we give estimates for the observability of the studied phenomena. 

[PACS numbers: 74.6O.Ge, 74.20.-z, 03.70. +k] 
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1 Introduction 



There has been a strong revival of interest, recently, in the physics of magnetic vortices in 
type II and high-temperature superconductors 0. The motion of these vortices subsequent 
to the applied supercurrent gives rise to a non-zero resistance, especially important in the 
case of the high-T c materials. Most research efforts have been devoted to phenomena relating 
to the nature of the mixed phase of a superconductor in some externally applied magnetic 
field and supercurrent. Issues connected with the pinning of the flux lines by defects have 
been widely studied. 

We 0, as well as Ao and Thouless 0] and Stephen M, have addressed the problem of 
the quantum dynamics of vortices in the absence of an external field but in the presence of 
an externally driven supercurrent, quantum dissipation and pinning. This leads to the decay 
of a supercurrent, or a residual low-temperature resistance in the superconductor. Whilst 
most of the dissipation seems to be ascribed to vortices tunneling in the sample from the 
edge 0, another interesting and novel possibility also explored by us in depth is that of a 
residual resistance arising from spontaneous vortex-antivortex pair creation in the bulk of, 
for example, a superconducting thin film. We have set up a powerful theory to study both 
bulk nucleation and tunneling from the edge of vortices spontaneously created either by the 
quantum fluctuations of the Magnus force field or by the current's gradient at the border of 
the sample. We find that, due to the high values of the critical currents in thin films, these 
dissipative effects can become particularly important in the two-dimensional (2D) geometry, 
where lateral fluctuations of the flux lines can be neglected. 

In this report we discuss our findings, comparing the rates of edge tunneling and bulk 
nucleation, to find that both effects can become important and observable under suitable 
conditions. Our theory is set up to study vortex tunneling and nucleation in the 2D geometry 
in the presence of quantum dissipation and of pinning potentials. The central result is that 
the tunneling/nucleation rate T has a strong exponential dependence on the number current 
density J 

Ttxv^ffJ- 1 exp{-77 e// 4/47rJ 2 } (1.1) 

Here r] e ff is an effective viscosity coefficient as renormalised by the magnetic-like part of the 
Magnus force, and £r is the barrier or nucleation energy for a single vortex as renormalized 
by screened Coulomb interactions and (fake) Landau-level corrections. This exponential 
dependence would make the tunneling/nucleation (folded, e.g., into the sample's resistance) 
observable in a rather narrow range of J-values. Thus, normally the superconductor is 
practically resistance-free. However, the high values of J that can be reached in the high-T c 
materials make the possibility of observing both tunneling and pair creation within reach 
for thin films. One particular feature that would uniquely relate the residual resistance to 
the phenomenon of spontaneous vortex-pair creation is the presence of oscillations in the 
J-dependence of T(J) in case a periodic pinning potential is artificially created in the film. 
These oscillations are in fact strictly connected to the pinning-lattice spacing d = 2ir/k of 
the periodic potential (we assume a square lattice), e.g. 

U(q(t)) = U j2[l-co S (kq a (t))] (1.2) 

0=1 

acting on the nucleating vortex-pairs described by a coordinate q. The problem of quantum 
dissipation for a particle (the 2D vortex) moving in a periodic potential has some interesting 
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features in its own right |5|, |6], [7|. It is characterised by a localization phase transition driven 
by dissipation; accordingly, two phases can occur depending on whether the dissipation 
coefficient || rj is greater (confined phase) or smaller (mobile phase) than a critical value 7] c = 
k 2 /2ir = 2ir/d 2 . This localization transition is described by a Kosterlitz-type renormalization 
group (RG) approach. We have implemented the RG approach for the evaluation of the 
dependence of the spontaneous nucleation rate of vortex-antivortex pairs on the external 
parameters for our own quantum dynamical system. A remnant of the dissipation-driven 
phase transition is observed and the pair production rate T can be derived in both phases 
by means of a frequency-space RG procedure leading to observable current-oscillations if 
i] > i] c . These result from a correction factor K(J) in Eq. ( p..l|) which is sensitive to the 
periodicity d of the pinning substate. By measuring these oscillations one can confirm the 
spontaneous pair-creation related nature of the residual resistance in a superconducting thin 
film, as well as microscopically or even mesoscopically probe the lattice structure of the 
pinning substrate. 

Below, we give a summary of our theoretical treatment and discuss the results obtained 
for the various cases relevant to supercurrent decay. Details of our work can be found in the 
literature 01. 



2 Theory of Spontaneous Vortex- Antivortex Pair Nu- 
cleation in the Presence of Quantum Dissipation 

One can show that a vortex moving in a superfluid experiences a force, the Magnus force, 
which is the equivalent of the electromagnetic Coulomb- Lorentz force. The classical equation 
of motion reads: 



mq = — V£7(q) + eE — eqxB — r?q (2.1) 

with m the (negligible) inertial mass of the vortex, carrying topological charge e = ±2ir, 
treated as a single, point-like particle of 2D coordinate q(t). U(q) is the phenomenological 
potential acting on the vortex and rj is a phenomenological friction coefficient taking dissi- 
pation into account. A supercurrent J gives rise to an electric-like field E = x J (a notation 
implying E • J = 0) superposed to a magnetic-like field B = zdp^ for a thin film orthogo- 
nal to the vector z and having thickness d with a superfluid component characterised by a 
3D number density The quantum- mechanical counterpart of Eq. ( |2.1| ) is constructed 
through the Feynman path-integral transposition in which the dissipation is treated quan- 
tistically through the formulation due to Caldeira and Leggett M. This approach views 
quantum dissipation as described by the linear coupling of the vortex coordinate to the 
coordinates of a bath of harmonic oscillators of prescribed dynamics. 

Quite generally, we deal with the quantum decay of the "vacuum" , represented by a static 
e.m.-like Magnus field, through its quantum fluctuations: in this case the vortex-antivortex 
pairs. The decay rate T is given by the formula (Af is a normalization factor) 

w fB = Im -AT(r) Vq(s)e-fo^ (2.2) 

2L U JO T Jq(0)=q(T)=q 

where L D , with D = 2 here, is the volume of the sample, and the Euclidean single-particle 
Lagrangian for closed trajectories in space-time is 
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(2.3) 



Here, the set {x^.} represents the Caldeira-Leggett bath of harmonic oscillators simulating 
quantum friction, with | Y^k mfc ^ fc ^(^— ^fc) = V u exp(— u/Q) (fl being some frequency cutoff) 
fixing the spectral distribution of frequencies (we restrict the present discussion to the ohmic 
case). Also, F^ u is the uniform Magnus field tensor and V(q) = 2£ U(q) is the relativistic 
counterpart of the classical potential, Eq. (|1.2j), thus having amplitude A = 2£ Uo (we have 
assumed Ao/£q<^.1). £q is the vortex nucleation energy. As one important remark, we assume 
here that the microscopic dissipation mechanism due to the normal material inside the vortex 
core, whether phonon- or electron-driven, can be always described phenomenologically by 
the Caldeira-Leggett type action. 

The Feynman path integral in Eq. ( |2.2| ) can now be evaluated at the leading singularity 
of the r-integral; in the dissipation-dominated regime (mi = m2 = m — > 0) we obtain the 
exact analytic expression for the bulk nucleation rate: 

L 2 32tt 2 J 1 ' J 

where T] e ff = \/rf + B 2 /rj and £q R = £q + f^Vv 2 + B 2 are renormalizations due to the 
magnetic-like field. Further renormalizations of £q due to vortex-vortex interactions are 
understood. According to this formula, for small current densities J the nucleation rate, 
hence the resistance [§], is negligible. However, depending on the material parameters, for 
higher J both the rate T and the film resistance can become important. Before discussing 
realistic estimates, we consider the startling effect of a periodic pinning potential on bulk 
nucleation. 



3 The Effect of a Periodic Pinning Potential: Current 
Oscillations 

Suppose a regular lattice of pinning sites is introduced artificially in the sample, e.g. in a 
square-lattice configuration. This can be achieved by means of a variety of etching or laser 
ablation or other techniques combined with modern MBE-type sample growth processes. 
We are then confronted with the motion (and then the quantum-mechanical nucleation) of 
particles in a dissipative and periodic-potential situation. This type of quantum diffusion has 
been studied at length [0. In the mobile phase, the RG analysis predicts a vanishing large- 
scale potential amplitude, thus no significant effect of the pinning potential can be expected 
on the nucleation rate, except for a further renormalization of £$. We look therefore for a 
situation corresponding to the confined phase (rj > r] c ) where the effective pinning potential 
amplitude grows up to an upper limit as nucleation proceedes in the bulk of the film p|. 

When the RG analysis is folded into the evaluation of the leading contribution to the 
Feynman path integral in Eq. (|2.2|) for the periodic pinning potential (|1.2j ), we end up with 
the following analytic result for the nucleation rate: 
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r = r K(j) (3.1; 

A V) = e(l + ,„)(l + ^)/„(^o(2H K : 

Here To is given by Eq. fl2.4j) , there is the further renormalization 8q R — > £q R + Aq for the 
nucleation energy, Jq{x) and Jo 0*0 are Bessel functions and the renormalised amplitude A 0R 
is a few times Aq itself 0. £n is a nucleation length, given in the lowest approximation by 
C-n = Sor/^ttJ. fj, is a dimensionless parameter of order unity. 

We then see that the competition between the nucleation length and the pinning- 
lattice spacing d = 2ir/k leads to oscillations, as the current J is changed, in the nucleation 
rate (hence in the resistance). This shows that in the nucleation process the dissociating 
vortex-antivortex pairs "feel" the periodicity of the pinning substrate. Now for the esti- 
mates of the expected values of T. So far, the computation referred to T = 0. In |2j we 
proposed to take into account finite temperatures by adding incoherently, in the above for- 
mulas, a term proportional to T (note ks = 1): J 2 ^J 2 + g^T, such that for J^O the 
standard Boltzmann exponential form in the vortex density, p v (x^/T, is recovered in the 
absence of pinning. This might be a reasonable qualitative treatment of thermal effects for 
temperatures less than ^8TrJ 2 /r]SoR- Furthermore, we propose, following Minnhagen ||, 
to account for vortex-antivortex Coulomb-like interactions by means of a current-dependent 
(and temperature- independent) activation energy: Sr(J) = £qr ln( J ma x/ J), with Sqr in- 
cluding all other renormalization effects. This expression entails an Sr infinite for J = 
(as is appropriate below the Kosterlitz-Thouless (KT) transition in the film) and changing 
sign, thus becoming unphysical, at Jmax- We take J ma x an order of magnitude higher than 
the single-crystal YBCO critical current, JL0| namely J^ x = 10 s A cm -2 , and take £r~80 
K for J = 10 7 A cm -2 , of the order of magnitude of the KT transition temperature M. 
The important friction coefficient is taken to be ^lO -2 A , from the Bardeen-Stephen 



formula WI\ applied to YBCO films. A most sensitive parameter is related to the amplitude 
of the pinning potential, e = Ao/£q R , which is entirely unknown. We have taken a negative 
e = —0.5, borrowing from classical nucleation the point of view that vortex production is 
actually aided by the presence of pinning centers. Finally, d = 50 A concludes our illustra- 
tive case which is clearly of the confined-phase type. As r\ > r) c , we have taken A 0R txA , for 
simplicity's sake, and find no significant dependence on the most uncertain of all parameters, 
Q. Its dimensions being the same as Sq, we have taken Q = £q R /c. Also, we take the film 
thickness s = 10 A as the typical interlayer spacing. Figure 1 then represents our estimate of 
the expected values of the nucleation rate in the above experimental conditions. Our claim is 
that both pair-production and current-oscillations should be observable with an appropriate 
technique (e.g. Tonomura's electron holography We point out that, with no controlled 

pinning present, the decay of a supercurrent in the absence of a magnetic field has been 



measured experimentally in BiSCCO films, see e.g. [T3|], though in this case edge tunneling 
(which we discuss below) also plays a role. 
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4 Theory of Vortex Nucleation from the Edge of the 
Sample 

Finally, we briefly describe our theoretical treatment for evaluating the tunneling rate from 
the edge (in the absence, however, of the periodic pinning potential). This is done by 
extending the treatment of Section 2 to the case where the vortices are already present in 
the edge strip of the sample. This requires that, in the formula for the rate, Eq. ( |2.2| ), the 
contribution from the antivortex paths is factored out and cancelled by the introduction of 
a suitable chemical potential /i*. Setting B = in order to simplify the treatment (a finite 
B results only in some renormalization), the formula for the tunneling rate becomes: 

r r°° dT r 

-wJm/ — $(T)/ V q (t)e- S ™ (4.1) 

2 JO 1 Jq(T)=q(0) 

where the non-relativistic (vortex-only) action reads 

-E.q+*7(q)J (4.2) 

and where the factor $(T) is fixed by a self-consistency procedure 0. The potential Z7(q) 
now contains the 2D Coulomb potential binding the vortices to the edge strip, a distance y 
away, and has the form U(q) = Ud(q) + K\n(l + y/a). Up is the usual Caldeira-Leggett 
potential and K the Coulomb interaction strength. At this point the usual leading contribu- 
tion to the T- and Feynman path-integral is extracted, however by means of a saddle-point 
approximation with inclusion of the Gaussian fluctuations M, which yields the tunneling 
rate per unit length (for an edge strip having width a): 

T_ = J e- So 

2L V a 2 ^/2Q^f{2TiJ-K/{a + y)) ^ ' ' 

Here, ip is a slowly- varying dimensionless factor of order unity, Q = ln[f2T/27r], and So, T 
and y are determined by the saddle-point conditions: 




_ 7r?7 _ 2 _ 77{inn(l + Kjl-nJa) + • • -} 2 

-2tt Jy + K ln(l + ^) = (4.4) 
a 

^-y - 2nJT + KT-^— = 
Q a + y 

This yields almost the same leading dependence (logarithmic corrections apart) r^e _ ^ Jo//J ^ 2 
on the external current that was also obtained for the homogeneous bulk-nucleation phe- 
nomenon. Indeed, the effective Coulomb energy i^ln(l + K/2irJa) can be interpreted as 
playing the role of the activation energy Sq renormalised by the vortex-antivortex interac- 
tions. As for the estimate of the edge tunneling rate, we use the same material parameters 
set up in Section 3, and in Figure 2 we report the tunneling rate R = T/L for different 
possible values of the Coulomb strength K. The conclusion is that both edge-tunneling 
(normally dominant for small J) and bulk-nucleation can become important and concomi- 
tant mechanisms for supercurrent decay in a superconducting thin film. 
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FIGURE CAPTIONS 



Figure 1. Vortex-antivortex production rate (in //m _2 s _1 , log 10 scale) in the presence 
(T) and absence (r ) of the periodic pinning potential, versus current density J em (in 
A cmr 2 x 10 7 ). K(J) is the current-dependent correction factor due to the periodic 
substrate. Here, the temperature is T = 2.5 K. 

Figure 2. Plot of the tunneling rate per unit length, R, as function of the supercurrent 
density, J em , for the values K = 30 K and K = 50 K. Dashed lines, T = K; full lines, 
T = 2.5 K. 
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